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I. INTRODUCTION 

The Kerr-Newman solution to the Einstein-Maxweh equations describes the gravitational field of a rotating and 
electromagnetically charged stationary black hole [1|. It generalizes both the static and charged Reissner- Nordstrom 
metric [2|, |3| as well as the rotating Kerr metric [J] . The latter is of very high importance not only for general relativity 
but also from an astrophysical point of view, as many black hole candidates were found in recent years, which are 
expected to rotate. Although it is not very likely that they also carry a net charge, some accretion scenarios were 
studied which may create such black holes 0, M, ■ 

On way to explore the gravitational field of a Kerr-Newman black hole is to consider the geodesic motion of (charged) 
test-particles in this space-time. Already shortly after the discovery of this solution many aspects of the geodesic 
motion were studied, among others, timelike equatorial and spherical orbits of uncharged particles [7| and the last 
stable orbit of charged particle s B (see also Sharp ^ and references within). Later the motion of charged particles 
was studied by Bicak et al [l3, lll| including, besides some general discussion for the radial motion, a discussion of a 
number of special cases like motion along the symmetry axis and circular motion of ultrarelativistic particles. Only 
recently, Kovaf et al [l3| found off-equatorial circular orbits of charged particles which are unstable outside the outer 
horizon, and Pugliese et al [ij] used equatorial circular orbits of neutral test-particles to distinguish between black 
holes and naked singularities. A comprehensive analysis of photon orbits in Kerr-Newman space-time was presented 
by Calvani and TuroUa [14] including the extended manifold with negative values of the radial coordinate and naked 
singularities. 

Analogously to the uncharged case, the geodesic equation in Kerr-Newman space-time can be separated by introduc- 
ing an additional constant of motion (besides the constants associated to the obvious symmetries of the space-time) , 
the Carter constant _15] . which ensures the integrability of the equations of motions. The resulting structure of the 
equations of motion is essentially the same as in Schwarzschild space-time, where they can be solved analytically in 
terms of elliptic functions as first demonstrated by Hagihara in 1931 [l6|. However, due to the remaining coupling of 
radial and colatitudinal equation, the generalization of his method to Kerr(-Newman) space-time was not straight- 
forward. This issue was solved by Mino [iTj ) by introducing a new time parameter, often called the Mino time, which 
completely decouples the equations of motion and enables a straightforward application of elliptic functions. This 
was already used to analytically solve the geodesic equation for bound timelike orbits in Kerr space-time by Fujita 
and Hikida [18] and for general timelike and lightlike orbits in Kerr space-time in [l9|. 

In this paper, we will discuss the geodesic motion of charged test-particles in Kerr-Newman black hole space- 
times. For the sake of completeness, we will include a magnetic charge of the black hole which was not done in 
the cited references but has interesting effects on the colatitudinal motion. After introducing the relevant notations 
and equations of motion in the next section, we proceed with a complete classification of timelike orbits of (charged) 
particles in Kerr-Newman space-time. This includes off-equatorial orbits, trajectories crossing the horizons, and orbits 
with negative values of the radial coordinate. In the fourth section, we will present analytical solutions in terms of 
elliptic functions dependent on the Mino time for all coordinates. The paper is closed by a summary and conclusion. 

II. GEODESICS IN KERR-NEWMAN SPACE-TIME 

A. Kerr-Newman space-time 

The Kerr-Newman spacetime is a stationary, axisymmetrical, and asymptotically flat solution of the Einstein- 
Maxwell equation 



G^. = -l(g''^F^^F,p ~ \g^,F^pF^^ 



where G^u is the Einstein tensor and F^^^ the electromagnetic tensor. Throughout, the units are chosen such that 
c = 1 for the speed of light and G = 1 for the gravitational constant. In Boyer-Lindquist coordinates the metric takes 
the form 

ds^ = P^dr^ + p'dO^ + ^ \{r^ + a^)dcj) - adtf - ^ \asm^e)dcl) - dt]^ (1) 

A p^ '- ■' p^ '■ ■' 



with 



p^{r,e)=r^ +a^cos'^9, (2) 

A(r) = r^ - 2Mr + a^ + Q^ + P^ , (3) 



and M > the mass, a the specific angular momentum, Q the electric, and P the magnetic charge of the gravitating 
source. (The existence of magnetic charges has not been proven yet but it will be considered for the sake of complete- 
ness.) We restrict ourselves here to the case that two horizons exist, given by the coordinate singularities A(r) = 0, 
r± = M zL \JM'^ — a? — Q"^ ~ P'^. The only genuine singularity is for p^ = 0, where r = and 9 = 7r/2 is fulfilled 
simultaneously. This means that a test particle approaching r = from above or below the equatorial plane does not 
terminate at r = as it would in Schwarzschild space-time but continues to negative values of r. For large negative 
values of r this can be interpreted as a "negative universe" , see {2(| . 

The Kerr- Newman metric reduces to the Kerr metric for Q = P — describing the exterior of rotating non charged 
black holes. It reduces to the Reissner-Nordstrom metric for a = which describes the exterior of a non rotating but 
charged black hole. In the case P = Q = a = Q the Kerr-Newman metric is reduced to the Schwarzschild metric. 

The electromagnetic potential is given by 

A = A^dx" = %{dt ~ asiT? ed(j)) + ^PcosOladt - (r^ + a^)d(j}) , (4) 

pZ pZ \ I 

from which the electromagnetic tensor can be calculated by F = ^{d,yA^ — d^A^)dx'^ A dx" . By the interchanges 
Q ^1' P, P ^> —Q the electromagnetic potential A of the dual electromagnetic tensor F can be obtained. 



B. Equations of motion 

The equations of motion for a test particle of normalized mass e, electric charge e, and magnetic charge h can be 
obtained by the Hamiltonian 



H = 7:9'"' {T^f. + P'K + hA^j){T:^ + eA^ + hA^) (5) 



where tt^ describe the generalized momenta. By introducing 

^ eQ + hP -- eP-hQ 



the Hamiltonian can be reduced to 



P^ . ^ , e=Ve^+h^ (6) 



H=-r''{TT^ + eA^){7T, + eA,), (7) 



where g and A are defined by ([T|) and ^ with P, Q, e replaced by P, Q, e. Therefore, the discussion of a test particle 
without magnetic charge is sufficient. In the following we omit the hats for brevity. 

We can obtain three constants of motion directly, since H does not depend on r, (p, or t. The first, e^ — —gufix'^x^^ is 
the normalization condition with e — \ for timelike and e = for lightlike trajectories. (The dot denotes differentiation 
with respect to an affine parameter r.) The second and third equation 

E = -TTt ^ -gtd - gt4,(t> + eAt , (8) 

L = n^ ^ g^ti + g,p,p(p - eA^ , (9) 

describe the conservation of energy E and angular momentum in z direction, respectively. A fourth constant of motion 
can be obtained by considering the Hamilton- Jacobi equation 

~drS = ^g'"'{df,S + eA,){d,S + eA^) . (10) 

With the ansatz S = ^er — Et + Lcj) + Si{r) + 5*2(0) it can be shown that the Hamilton- Jacobi equation indeed 
separates with the Carter constant K as separation constant, see |15| . 



With these four constants the equations of motion become 

(f.) ^K-ea^cosH-^^^:e{e), (11) 

\djj sin^e ^ ' ^ ' 



(^^) =7^2(r-)-(er-2 + if)A(r-)=:i?(r-), (12) 



dcj) _ aR{f) 7(9) 
dj ~ A(f) sin^ e ' 
dt (f2+a2)7^(f) 



(13) 
aTie) , (14) 



dj A(f) 

where 

n{f) ^ {f^ + a^)E - aL - eQf , (15) 

Tie) ^dEsin^ 9 ~L + eP cose. (16) 

All quantities with a bar are normalized to M, i.e. a; = xAI for a; = r,a,L,Q,P as well as /iT = KAP and A(f) = 
f^ — 2f + a^ + Q^ + P^- Here 7 is the normalized Mino time [13] given by d-y = Alp^^dr with the eigentime t. 

III. CLASSIFICATION OF MOTION 

In this section we will classify the types of orbits in terms of colatitude and radial motion. We will analyze which 
orbit configuration, i.e. which set of orbit types, may appear for given parameters a, Q, P, E, L, K, e and which region 
in parameter space a given orbit configurations occupies. Here we assume that e = 1, that is, we restrict ourselves to 
test particles with mass, but the discussion for light may be done analogously. The whole analysis will be based on 
the conditions Q{9) > and R{f) > which are necessary for geodesic motion. 

For both colatitudinal and radial motion, we will first give some general properties as symmetries and notation of 
orbit types. We then proceed with the determination of possible orbit configurations, i.e. sets of orbit types which 
are possible for given parameters (for more than one possible orbit type the actual orbit is determined by initial 
conditions). Each orbit configuration covers a particular region in the parameter space. Finally, we will analyze how 
these regions look like and determine their boundaries in parameter space. 

A. Colatitudinal motion 

The coordinate 9 may only take a specific value 6*0 G [0, tt] if 0(6*0) > is valid. We will analyze in the following 
whether this is fulfilled for a given parameter set. 

1. General properties 

First, we notice that Q does not depend on Q and that e and P only appear as the product eP. The function Q 
has the following symmetries: 

• A change of sign of eP has the same effect as reflecting 9 at the equatorial plane: Q\_^p{9) — Q\^p{tt ~ 9). In 
particular is Q symmetric with respect to the equatorial plane if eP = 0: Q\eP=o{9) — Q\i,p=oi'^ ~ ^)- 

• A simultaneous change of sign of Z and i? result in a reflection at the equatorial plane: 0|_£ -~e{9) = '^\l si'^^^) 
ande£^o=B(0) = e£^o^^(^-0). 

Therefore, we assume w.l.o.g. eP > and E > in the following. The condition 9(6') > also shows that A" > is 
a necessary condition for geodesic motion as all other (positive) terms are subtracted. 

The Carter constant also encodes some geometrical information if considered in its alternative form C — K ~ {aE — 
LY . Because of 8(7r/2) = K — [aE — L)^ = C a particle may only cross or stay in the equatorial plane if (7 > 0. 
For equatorial orbits even (7 = is necessary as 3^ (f ) =0 needs to be fulfilled. For K = geodesic motion is only 
possible if a = L = or in the equatorial plane with dE = L. 

The function 9(6*) constains a term which diverges for 9 — 0,n given by ~ l^in^'e"'^ • This fact suggests to 
distinguish between two cases: 



• L 7^ zLeP: In this case Q{9) -^ — oo for 6* ^ 0,7r, that is, the north or south pole will never be reached. 

• L — ±eP: In this case 8(6') —i'K — a^ for ^- and L = eP as well as for 6* — ^ tt and L — —eP. Therefore, 
a particle with L = eP may reach the north pole and a particle with L = —eP the south pole if in addition 
K > a^. For the subcase L = = eP both north and south pole may be reached if ^ > a?. 

If the parameters are such that the poles 6* = 0, tt can not be reached it is convenient to consider ly = cos 9 instead of 
9. In terms of v the differential equation for colatitudinal motion reads 

^)'.^6,.'^=:e.(.), (17) 

where bo ^ K - {L - aEf = C, h = 2eP{L - aE), b2 = -K - a^ - e^P^ + 2a^E^ - 2aEL, 63 = 2aEeP, and 
64 = 5^(1-^^). 

For a given set of parameters of the space-time and the particle different types of orbits may be possible. We call 
an orbit 

• northern or N, if it stays in the northern hemisphere 9 < n/2, 

• normal or E, if it crosses or stays in the equatorial plane 9 = 7r/2, 

• southern or S, if it stays in the southern hemisphere 9 > 7r/2. 

Equatorial orbits with 9 = tt/2 are a special case of normal orbits. In addition to the above notions we add an index 
N if the north pole 9 = and S if the south pole 9 = tt is reached, for example Nn for a northern orbit reaching the 
north pole. 

2. Orbit configurations 

Let us now analyze which orbit configurations, i.e. which sets of the above introduced orbit types, are possible 
for given parameters. We use the necessary condition for colatitudinal motion Q{9) > for this, which implies to 
analyze the occurrence of real zeros of 9 in [0, tt] and the behavior of Q at the boundaries = 0, tt giving the sign of 
Q between its zeros. (The orbit a test particle with the given parameters actually follows in a space-time with the 
given parameters depends on the initial values.) 

(A) Case L 7^ zLeP Here Q{9) = —00 at = 0, tt which implies that Q^, has an even number of zeros in (—1, 1) 
(counted with multiplicity). If Q^, has no real zeros there, no colatitudinal motion is possible, which gives a restriction 
to the permitted sets of parameters for geodesic motion. In the case of 2 real zeros there is a single orbit of type N, 
E, or S, which is stable at a constant 9 if the 2 zeros coincide. If all 4 zeros of &i, lie in (—1, 1) all combinations of 
two orbit types except EE are possible. For two or more coinciding zeros this point is stable if it is a maximum of 
&1, and unstable otherwise. 

(B) Case L = ±eP Here we have to consider four subcases: (Bl) K < a?: The same orbit types as in (A) are 
possible. (B2) K > a^ , L ^ 0: In this case G has different signs at = 0, tt which implies that O has an odd number 
of real zeros in (0, tt). For 1 real zero there is one orbit of type Nn or En for L = eP and one of type Ss or £^5 for 
L = —eP. If Q has 3 real zeros in (0, tt) there is one additional orbit not reaching a pole. (B3) K > a^, L = = eP: 
Here > at = 0, tt, i.e. Q has an even number of zeros in (0, tt). Also, for eP = the function Q is symmetric with 
respect to the equatorial plane. For no real zeros there is one orbit of type -Ejvs which reaches both poles, and for 2 
real zeros an Nn and an Ss orbit. More zeros in [0,7r] are not possible. (B4) K = a^: Here an orbit with constant 
9 = {9 = tt) is possible for L = eP [L = ~eP). For Z = no other than the two constant orbits are possible but 
for Z 7^ it is 6 — ?► —00 at the other boundary. In the latter case, if the orbit is stable, there may be one additional 
orbit of type E, N, or S. 

For an overview of this different orbit configuration see table HI 

3. Regions of orbit configurations in parameter space 

It is now of interest for which sets of parameters a given orbit configuration changes. As > is necessary for 
geodesic motion, this happens if the behavior of Q at the boundaries changes, which means a switch from one of the 
above cases (A), (Bl), . . . , (B4) to another, or if the number of real zeros of Q changes. The latter occurs at that 
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TABLE I: Overview of different orbit configurations for colatitudinal motion. The vertical bar of the second column 

denotes 6 = tt/2 and the thick lines G > 0, i.e. regions where a motion is possible. Dots represent single zeros and 

circles double zeros. If zeros merge, the resulting orbits are stable if a line is reduced to a point and unstable if lines 

merge. The configurations {B2) and {B4) with L — —eP are obtained by a reflection at the equatorial plane. 



parameters for which Q has multiple zeros. With these two conditions the different regions of orbit configurations in 
parameter space can be completely determined. The first condition was already analyzed above. 

For 9 € (0, tt) the function 9 has the same zeros as the polynomial 9^ and we may use 9,^ instead of 9 for all 
orbits not reaching — 0, tt. If L = ±eP then vq = ±1 is a zero of 9^ but does not correspond to a turning point of 
the colatitudinal motion. If in addition K = a? then vq ~ ±1 is a double zero of 9,^ and = 0, tt a simple zero of 9, 
which does correspond to a turning point of 6. Keeping this in mind we will also use 9^ for these cases but discuss 
the occurence of multiple zeros at = 0, tt separatly using 9. 

The condition for a double zero vq is ^f (t'o) = = 9^(t'o). This can be read as 2 conditions on 2 of the 5 
parameters E, L, K, eP, and a. Solving these two conditions for E and L dependent on the position of the double 
zero vq and the other parameters yields 
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(18) 
(19) 



The expressions for E and L diverge at i^o == for eP ^ 0, at vq = \J K /a?, and at i^g = ±1 for E, what suggests to 
consider the 2 conditions for double zeros directly for these points: 

Equatorial orbits For J^o = the 2 conditions imply that either K = and L — aE oi K — [dE — L)^ and eP — 
are necessary and sufficient for the existence of equatorial orbits. The asymptotic behavior of E and L &t vq = also 
displays these conditions, 
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eP 
2diyf) 
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(20) 
(21) 



In the case K = 0, L — aE the equatorial orbit is the only possible geodesic orbit (for G ^ 0) and, thus, stable, 
whereas for K = {aE — L)^, eP = the sign of A := d'^{E^ — 1) — Z^ has to be considered: The orbit is stable ii A < 
(for 8^0) and unstable if ^ > 0. In the case A — the zero 9 — tt/2 is even fourfold. For K = a^, eP — — L, and 
E'^ = 1 the function 8 is identical to zero. 

Orbits with 9 = Q,tt Geodesic motion along the axis = 0, tt is possible only if 6* = 0, tt is a double zero of 8. The 
2 conditions on double zeros show that L — eP and K — a? are necessary and sufficient for 9 = 0, and L = —eP, 
K = a? for 9 = TT. This can also be seen by considering the asymptotic behaviour of Li,2 and Ei,2 as vq approaches 
±1: it is given by lim^Q_j.i L1.2 — > eP and limj,(,_j._i Li_2 = —eP whereas £^1.2 diverges for K 7^ a?. 

Let us discuss the stability of the orbits 9 = 0,7t: The orbits are unstable if a^ — {aE — L/2)'^ > and stable if 
a^ — {aE — L/2)^ < 0. For a^ — {aE — L/2)^ = the poles 9 = 0,Tr are fourfold zeros and the orbit 9 = is stable 
if _E = ^ ~ 1' ^nd the orbit 9 = niiE = ^ + l. In the special case of L = = eP the two orbits are unstable if 
£;2 < 1 and stable if E' ^ > 1. For L = = eP, K = a'^, a nd E'^ = 1 again 8 = 0. 

Orbits with 9 = ±\/K/a'^ The singularity i/q = ±y^K/a'^ is located in (—1, 1) and is not equal to zero only if 
< K < a^. Assuming this, orbits with constant cos 6' = ±y/K/a'^ can exist only if aL — {a^ — K)E = ±ePvK is 
fulfilled. This can be infered from the asymptotes of Z1.2 in terms of £'1.2 around i^o = iy/K/a^, 

a^ -K ^ eP^K 



^1,2 = ^£i,2±^^ + V^TA/r^ . (22) 

a a \ \ y a'' j 

The expressions P^ and P^ depend linearly on eP and a rescaling of parameters {L/a, K/a"^, eP/a) removes 
the rotation parameter a completely from the equations. Only the dependence on K is not obvious. If we solve the 2 
conditions ^f(^o) = = G)u{^o) for E and K instead of L this yields 
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2 eP{iy^ + 1) - 2iyoL 



24{vl - 1) 



{eP{v'^ + 1) - 2voL) ± J{eP{yl + 1) - 2LvoY + H^^i^l - 1 



(24) 



Note that Ei and Ki behave regular at i^o = and aEi = L + O{i'o), Ki = 0{v^) whereas Ei and K^ are finite 
at 1^0 = only if eP = 0. Then K — >■ {aE — Z)^ as expected from the analysis of equatorial orbits above. Also, for 
vq — > ±1 the conditions L = ±eP, K ^ a} are recovered. 

Triple zeros are also of interest as they correspond to parameters where a stable orbit with constant 9 becomes 
unstable and vice versa. Therefore, will will study them here. The 3 conditions for such points are = '^^■2 {vq) = 
^f (j^o) = '^v{vq) which we read as 3 conditions on E, Z, and eP yielding 

^ _ ^ 1 2vl-a^ + &vlK-a^ - (3^.0^5^ + g)(g + a^) 
''' 2 V(l - i^'^){K - a2^o')(l - ^.i){K - -a^ulYa ' 



1 v/(l - vl){K - a?vl){K - a^){(Pvl - i v^a' + ivlK - K) 

2 (^2 _ 1)2(^.^2^2)2 






eP„ = ±^ . (27) 

V(l - ^2)(A- - aV„){i.l - 1)(A- - oVg) 

In particular, triple zeros are also double zeros and the asymptotic behavior of Z as a function of E and K at the 
singularities vq = ±^JK /a^ is determined by (|22|) . The points i^o = ±1 were studied above but we note here that 
= 0; vr is a double zero of 8 if i^ = ±1 is a triple zeros of 81,. 

Besides the equatorial orbits and orbits with constant = 0,7r discussed above, fourfold or even higher order zeros 
are only possible for parameters corresponding to 8 = 0, which are a = eP = L = K = with arbitrary E ot E = ±1, 
K = 52, and Z = eP = 0. 

The results of this section are visualized in figures [T] to |4l As the five dimensional parameter space can not be 
completely pictured we have to fix at least two parameters. For this we choose a, which can be removed from equations 
(UHl), (ITOl) . (P5)) . and (PH) by a rescaling of parameters, and eP, which enters only linearly in P^ an ([TO|. As a three 
dimensional plot of E, Z, and K is still confusing we present two dimensional plots of Z over E and K over E. Special 
cases to be discussed are then eP = 0, Z = ±eP, K = a? , and ^ = 0. For the latter geodesic motion is only possible 
on a stable orbit of constant 9 (if not 8 = 0), more precise on the equator for Z = aE or on ePcos9 = Z for a = 0. 
In each plot, we use the following conventions: 
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FIG. 1: Orbit configurations for the colatitudinal motion with a — 0.5, K > a^, and eP > 0. For a general 

description see the text. The orbit configurations on the sohd lines contain an orbit with constant 9. They are 

unstable if marked by the red solid line starting at the dot and approaching the dashed line, and stable 

otherwise. The dash dotted line a.t L — eP (L — ~eP) denotes an orbit crossing the north (south) pole and 

corresponds to the orbit configurations {B2). The other regions correspond to the configurations {A). Small plots on 

the top are enlarged details of the lower plot. 

• Solid lines indicate double zeros of Q which correspond to stable or unstable orbits of constant 0. We use red 
lines for orbits with constant 9 < ^, blue for constant ^ > §, and green for equatorial orbits. 



• Dashed lines denote orbits with turning points at the equatorial plane and are given by Z = aE ± v K . They 
also mark the transition from N or S to E orbits and are asymptotically approached by solid lines corresponding 
to orbits with constant 9 near the equatorial plane. 

• Dash dotted lines correspond to orbits which cross a pole. They only appear in L over E plots for K > a^ 
and are located at L = ±eP. In this case they are asymptotes to the solid lines corresponding to orbits with 
constant 9 near a pole. 

• Dotted lines are asymptotes to solid lines. They do NOT separate different orbit configurations. For L over E 
plots they only appear for K < a? and are approached for fg -^ -iivK /a. If K is plotted versus E they are 
approached for vq — > ±1. 

• Single dots mark triple zeros which separate stable from unstable orbits. 

• The labels iV, E, and S indicate the orbit configurations summarized in table HI Solid, dashed, and dash dotted 
lines indicate transitions from one orbit configuration to another. Here dash dotted lines are special in the sense 
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FIG. 2: Orbit configurations for the colatitudinal motion with a = 0.5, K < a^, and eP — 0.3. For a general 

description see the text. The orange solid lines in (a) correspond to orbits of constant 6* = 0,7r, which are unstable 

for energies less than the one marked by the dots. These orbits correspond to the configurations (-84), other regions 

(in both plots) to the configurations (A) or (Bl). In both plots, the orbits labeled by the red solid line starting at 

the dot and approaching the dashed line are unstable. All other orbits on solid lines are stable. Small plots on the 

top are enlarged details of the lower plot. 



that on both sides there are always the same orbit configurations but only the line itself corresponds to another 
configuration. Regions where a colatitudinal motion is forbidden are marked with 'none'. 



Note that we plot usually only positive values of E as negative values can be obtained by E 
reflection at the equatorial plane (e.g. an A^ orbit becomes an S orbit). 



^E, L 



-L and a 



B. Radial motion 

The discussion of the radial motion will be analogous to that of the 0-motion. An orbit can only take a specific 
value in [— oo, cxd] if the radicand of the right hand side of eq. (J12p is larger than or equal to zero at that point, i.e. if 

R{f) = {{f"^ + a^)E - al - eQff ~ {ef'^ + K){f^ - 2f + a^ + Q^ + P'^) > 0. 



We will analyze for which values of the parameters E, L, K, e, P and Q this inequality is satisfied. 
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-2-lr- 





(a) K = 0.5 



(b) K = 0.2 



FIG. 3: Orbit configurations for the colatitudinal motion with a — 0.5 and eP — 0. For a general description see the 

text. The green sohd hnes correspond to equatorial orbits which are unstable from dot to larger L and else stable. 

The black solid line marks two stable orbits of constant 7^ § which are symmetric with respect to the equatorial 

plane. The dash dotted line in (a) marks the case (-B3) with an Ens orbit between the green solid lines and an 

NpfSs configuration else. 



General properties 



The three parameters e, Q, and P appear only in the two combinations eQ =: Q and Q^ + P^ =: V^ which may be 
considered instead. The function R has the following symmetries: 

• It depends quadratically on P (and V): -R|_p — R\p {R\-v = R\v)- 

• A simultaneous change of sign of Q and e results in the same motion: R\^q _e = R\q e- 



• Also, a simultaneous change of sign of L, E and e(or Q) results in the same motion: R\-e.-l - 



R\ 



E,L,e 



Thus it suffices to consider P > 0, > 0, and e > (or, equivalently, V > and Q > 0). Note that ^ > was a 
necessary condition for the colatitudinal motion to be possible at all and, therefore, this condition remains valid. 

The sign oi C = K — {aE — L)^ again encodes some geometrical information. At f = the polynomial -R(O) — 
-(Q^ + P'^){L - clEY - (a^ + Q^ + P'^)C can only be positive if (7 < 0. Since C" > needs to be satisfied for an 
orbit to reach the equatorial plane this implies that C = is necessary for an orbit to hit the ring singularity. For 
q2 ^ p2 ^ Q g^jgp j^ _ g^ = is neccessary to hit the singularity and, thus, also K = Q. 

The zeros of the parabola A(f ) are the horizons f± — 1± \/\ — a^ + Q-^ + P^ with A(f) < in between. Therefore, 
the polynomial R(f) is always postive for f g [r_,f+]. This implies that there can not be any turning points or 
spherical orbits between the horizons. There is a turning point at a horizon if 7?.(f-|-) = (see eq. (1151) ). 

For a given set of parameters of the space-time and the particle different types of orbits may be possible, for which 
we use the terminologies 

• transit or T, if f starts at ±cx) and ends at Too, 

• fiyby or P, if f starts and ends at +00 or — cxd. 



bound or B, if r remains in a finite interval [r 



mm; ' maxj 



We add an index +, 0, or — to a flyby or bound orbit if it stays at f > 0, crosses f = 0, or stays at f < 0. Also, a 
superscript * will be added if the orbit crosses the horizons, i.e. contains the interval [f_,f+]. For example, the orbit 
F'^ comes from infinity, crosses the two horizons, turns at some < f < f_, and goes back to infinity. If more than 
one orbit is possible for a given set of parameters, the actual orbit of the test particle is determined by the initial 
conditions. 
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-3.5 -3 -2.5 -2 
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(a) L = 0.5 



(b) L = eP 



FIG. 4: Orbit configurations for the colatitudinal motion with a = 0.5 and eP = 0.3. For a general description see 

the text. For both plots, solid lines with dashed asymptotes correspond to unstable and all other to stable orbits. In 

(b) the regions above the solid orange line at K — cP correspond to the configurations {B2) and below to (Bl). The 

line itself corresponds to (54) and covers all possible configurations (leftmost corresponds to topmost in table |T]) 

with the exception of the N^^N configuration which is only possible for large values of eP. 



2. Orbit configurations 



We will now analyze which sets of the above introduced orbit types are possible for given parameters. Geodesic 
motion is only possible in regions of R{f) > and, therefore, the possible orbit configurations arc fully determined 
by the number of real zeros of the polynomial R and its sign at zLoo. The latter is determined by the sign oi E^ ~ 1 
(or the lower order coefficients ii E^ ~ 1, what we will consider separately). This suggest to introduce the following 
classes of orbit configurations. 

(I) Case E^ > 1 Here R{f) — > oo for f — >■ ±cx) and R may have none, two, or four real zeros. For no zeros there 
is a transit orbit, for two zeros there are two fiyby orbits, and for four zeros there are two flyby and a bound orbit. 

(II) Case E^ < I For such energies a test particle can not reach ±oo, R{r) — > — oo for f — > icxi. As R{f) > 
between the horizons it has at least 2 real zeros and there is always one bound orbit crossing the horizons. If R has 
four real zeros there is an additional bound orbit. 

(III) Case E^ ^ 1 Here the behaviour of R at infinity depends on the sign of 1 — eQ. For 1 — eQ > it is 
R{f) — ?> ±oo for r — > ±oo and the other way around for 1 — eQ < 0. In both cases R has one or three real zeros. For 
one real zero there is a flyby orbit which crosses the horizons and for three real zeros there is a flyby and a bound 
orbit. The flyby orbit reaches +C!0 for 1 — eQ > and — oo for 1 — eQ < 0. If also 1 — eQ = one has to consider the 
sign of the second order coefficient and so on. 

For an overview of this different orbit configuration for the radial motion see table [Til 
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real zeros 


range of r 


types of orbits 







T 




2 
2 
2 
2 


1 , , « - 


F-,f; 

F-,Fi; 










4 
4 
4 
4 
4 
4 
4 
4 
4 




Fo, B+, F+ 
Fo, Bl, F+ 
F-,BX,F+ 
F-,B5,F+ 
Fo, B+, Fl 
F-, B+, F; 
F-,Bo,Fl 
F_, B_, Fl 
F-,B-,FS 





















(I) E'^ > 1 



real zeros 


range of r 


types of orbits 


1 








Fo 








3 
3 
3 
3 








F-, BS 
F-,Bl 
Fo, Bl 
F5, B+ 



























real zeros 


range of r 


types of orbits 


2 
2 


I -. , -. 


Bo 






4 
4 

4 
4 
4 
4 




B+, B+ 
B'o, B+ 
B+, BX 
Bo, B'X 
B-, BX 

B-, b:, 
















(II) £2 < 1 





real zeros 


range of r 


types of orbits 


1 
1 




FX 
FX 






3 
3 
3 
3 
3 
3 




BX,F+ 
Btu F+ 
B+,FX 
Bo, FX 
B-,FX 
B-,FX 




1 A •- - ' 











(III) E'^ = 1, eQ < 1 



(III) E^ 



1, eQ > 1 



TABLE II: Overview of different orbit configurations for radial motion. The vertical bars of the second column mark 
f = 0, f ~ f^, and f = f+ (from left to right). The dots represent the real zeros of R (turning points) and the thick 
lines R> 0, i.e. regions where a motion is possible. If zeros merge, the resulting orbits are stable if a line is reduced 

to a point and unstable if lines merge. 

3. Regions of orbit configurations in parameter space 

After considering the possible orbit configurations for the radial motion we will now find the sets of parameters for 
which a given orbit configuration changes to another. The orbit configurations are fully determined by the signs of 
i?(±oo) as categorized above and the number of real zeros of R, which changes if two zeros of R merge. The latter 
occurs if the two conditions on double zeros R{fo) = and ^(^^o) = are fulfilled. Read as two conditions on E and 
L this implies 



^ _ eQ ^ 1 VA(r-o)(fg + i^)(foA(fo) + (fp -l)iK + F^ 
''" 2fo 2 A(fo)(r2+i^)fo 



Li, 



eQia' - f2) 1 VA(fo)(f2 + if )(foA(fo) + (fo - l)iK + fg))^ 



± 



2afo 2 A(fo)(f2 + K){foA{fo) + (fo - 1){K + f^))foa 

X (foA(fo)(a2 ^K) + {K + fl){fl - f„{P^ + Q2) _ ^2)) ^ 



(28) 



(29) 



These expressions diverge at tq = and at the horizons rg = r±. We consider these points separately, see below. In 
addition, L diverges also at fo — ±oo, Z1.2 — > ^^en{ro)-eQ -^ ^ 0(1) for fo -^ ±00. This implies that for eQ > 1 
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other orbit configurations than for eQ < 1 may appear. At the limits fo —> ±cx3 the energy E remains finite, 
Ei_2 — i> ±sign(fo) there. 

Orbits at f — As the ring singularity is located at f = 0, 9 = ^ only those orbits which are not equatorial do 
not terminate at f = 0. From the discussion of the colatitudinal motion equatorial orbits occur for K — 0, aE = L or 
K = {aE — L)^, eP = 0. If we exclude these parameters f = is a multiple zero if and only if eQ y^ 0, K — e^Q^A{0), 

and aE — L — Si£_LJ., This can also be infered from the asymptotic behaviour of Ei^2 and Zi_2 at tq — 0, which is 
given by 



""^'^ ^ — ^aWo — ^ ^^om + ^ °^ ' ^ ^ 

Li,2 = ai?i.2 ± ^^-^5^ + O(ro) . (31) 

a 

Here the equation for E2 implies two facts: First, at .R" = e^Q^A(O) the regions of orbit configurations essentially 
change because the sign of the r^^ term changes and £'2 switches between ±00. Second, for K = e^Q^A{0) the 
expressions for E2, L2, and eQ are the same as p5l) - ((37)) with fo = 0, which correspond to triple zeros. 

Orbits at the horizons It was noted above that a horizon is a turning point of the radial motion if 'JZ{f±) — 0. This 
implies that the horizons can not be multiple zeros because of ^(f±) = — 2(f± — l){f^ + K) 7^ for TZ{r±) = 0. (This 
is valid for massive test particles only. For light rays a horizon is a multiple zero if A" = and aL = E{f'j. + a^) — eQf± 
but the corresponding orbit is always unstable.) The asymptotic behaviour of Ei 2 and Li 2 near the horizons is given 

by 



72 \/^ 


-V' 


% + K 


4 f+ 


Vo - 


-T+ 


V2V^ 


-IWf 


-■i+K 



lim Ei;2^±-, , " +0{1), (32) 



lim Ei,2^T-, _ ,_ _ +0(1), (33) 

\imLi2^—i2f±-P^-Q^)+Oil). (34) 

ro—>r± a 

Now let us turn back to the general case. The expressions ([^51) and (P5)) depend linearly on Q = eQ and the 
parameter a can be removed by a rescaling of the other parameters (4, -S, S^, i|-) and the radial coordinate (?■). 
(Note that also distances measured in units of M have to be rescaled). The dependence on K and V^ = Q^ + P^ is less 
obvious. This can be studied by solving (j28p and (j29p for these parameters but the expressions are quite cumbersome 
and we do not give them here. However, in Figure IH] orbit configurations for varying K are shown and all possible 
orbit types already appear in Figures [5] to [3 

Let us now analyze where triple zeros occur as they mark transitions from stable to unstable orbits. Solving the 
three conditions ^{fo) = ^(^0) = -R(^o) = for E, L, and Q ^ eQ yields 

^ _ ^ 1 KA^(fo) + (fg + K){3fl 2r-o + K)A(fo) (fg + K)^fo - if 
''" 2 ifl+K)Aifo)V{rl + K)A{fo) 

L,,2 = T^ X^i°^,Syr [^(3-1 - «' + 2i?)A2(fo) + {fl + K){ft - Kifl - 2f, + a') ~ aHSfl - 2fo))Aifo) 

+ ifl + K)^{-a^-fl)if,-lf], (36) 

Q _ fgA^(r-o) - (fg + jr)(2r-g - fg + J?)A(fo) + (fg + K^foifo - 1)' 

ifl + K)Aifo)Vir'o+K)Aifo) 

Here again, the horizons f± are singularities and in between triple zeros are not possible. Fourfold zeros may only 
occur for f < 0, given by E,,2 = ± ^ Y'\ , L,,2 = ±^^^20^^ Q ^ ±^=f==, and K = ^^4^' or at the 

ring singularity HE — ± ,} , A' = 0, 2=4, and Z — aE. 

tot. J ^A(O) -^ 

The regions of different orbit configurations are visualized in Figures [5][8l As the parameter space is six dimensional, 
we have to fix at least three parameters for plotting. We always choose to fix a and Q = eQ for the reasons outlined 
above. Here Q = ±1 and Q = ±y^K/A{0) will separate quite different plot structures: at Q = ±1 the behaviour of 
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4- F 




0.8 0.9 1.0 1.1 1.2 1.3 1.4 2 3 4 5 6 7 8 




(a) eQ = 0.3, eQ < 1< J -^ 



(0) 




FB^F^ 



B,B 



+ + I 



0.5J — ^ 



0.95 1.00 1.05 0.96 0.98 1.00 1.02 1.04 

-1.7i, 




-1.1 -1.0 -0.9 -0.8 -0.7 -0.6 




(b) eQ = 0.9, eQ < 1< ^^ 

FIG. 5: Orbit configurations for the radial motion with a = 0.6, K = 1, and V^ — 0.4 for case (i). For a general 

description of colours and linestyles see the text. The spherical orbits marked by the dark blue (green) solid lines 

starting at the dots and approaching the light blue dashed (the black dash dotted) lines are stable. In (b) also the 

orbits on the red solid line between the two dots are stable. All other spherical orbits are unstable. Note that for the 

small eQ of (a) regions of orbit configurations are only slightly deformed by the transformation E — > —E, L — > — L 

but remain unchanged otherwise. Small plots on the top are enlarged details of the lower plot. 



Li^2 at infinity changes and at Q = ±y^K/A{0) the behaviour of i? at r = 0. We will restrict here to positive values 
of Q but allow all values of L and E. The values eQ < are recovered by {E, L) -^ {—E, — L) as noted in the section 

on general properties. Therefore, we distinguish between four different regions: (i) Q < minjl, W . , , }, (ii) < / . , . < 
Q < 1, (Hi) 1 < Q < J -Sq), and (iv) max{l, J -S^} < Q- Note that for eC, 



the polynomial R is unchanged 
{—E, —L) and, therefore, regions of orbit configurations differ only slightly from this 

^(A(0) 



by the transformation {E^ L) 

symmetry for small Q in region (i). However, if Q is larger than (j37p with tq 



K - y/ {/^{Q) - Kf + AK) 
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IM. 



F-B+F+ 




0.96 0.98 1.00 1.02 1.04 



-1.5 




^jy < eg < 1. 



FIG. 6: Orbit configurations for the radial motion with a = 0.6, K = 0.5, V^ = 0.4, and eQ = 0.9 

For a general description of colours and linestyles see the text. The spherical orbits on solid lines starting at a dot 

and approaching the black dash dotted lines or the light blue dashed lines are stable (two green lines, a dark blue 

line, and a red line). All other spherical orbits are unstable. In the detailed plot on the upper right the solid blue 

line approaches the light blue dashed line so close from below that they are hard to distinguish; the regions 

indicated there are meant to be between them. The same holds for the lower right plot with the red solid line 

approaching from above. 



(local minimum) the structure changes as an additional pair of triple zeros appears, see Figure [5bl 

For better comparison, we always plot L over E and indicate the dependence on K and V — P^ + Q"^ hy slowly 
varying them in a plot series. In each plot, we use the following conventions: 

• Solid lines indicate double zeros of R which correspond to stable or unstable spherical orbits of constant f. We 
use red lines for orbits with constant f < 0, blue for constant < f < f_, and green for constant f > f+. 



• Dashed lines denote orbits with turning points at f = and are given by L± = aE ± ^ — ^-^. They mark 
transitions between orbits with different indices. Red or blue solid lines corresponding to orbits with constant 
r near r — asymptotically approach L±. 

• The dash dotted line marks E^ = 1. Between E = ±1 no orbit can reach infinity. In addition, £■ = ±1 is 
asymptotically approached by red or green solid lines with r — >■ ±oo 

• Dotted lines are asymptotes to solid lines for f —>■ f±. They do NOT separate different orbit configurations. 

• Single dots mark triple zeros which separate stable from unstable orbits. 

• The labels T, F, and B indicated the orbit configurations summarized in table HT] Solid, dashed, and dash 
dotted lines indicate transitions from one orbit configuration to another. 
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2 4 6 8 




FIG. 7: Orbit configurations for the radial motion with a = 0.6, V^ = 0.4, K = 1, and eQ — 1.1, 1 < eQ < \/ xf^:- 

For a general description of colours and linestyles see the text. The spherical orbits on solid lines starting at a dot 

and approaching the black dash dotted lines or the light blue dashed lines are stable (two dark blue lines, a green 

line, and a red line). All other spherical orbits are unstable. In the detailed plot on the upper right one of the solid 

red lines approaches the light blue dashed line so close from below that they are hard to distinguish; the regions 

indicated there are meant to be between them. 



IV. ANALYTICAL SOLUTIONS 

We will now solve the equations of motion (fTTI) - ([T^ with the initial conditions 

^(70) = 00 , ^(70) = fo , 0(70) = 00 , ^(70) = to , (38) 

In addition the initial direction, i.e. the sign of ^(70) for x — 9,f, has to be specified. We denote this by a^ = 
Bgn(^(7o)). 



A. ^-motion 

The equation of motion (|11|) needs to be solved. We first concentrate on the case in which the poles are not reached, 
i.e. 0(7) € (0,7r). Then it is convenient to substitute u — cos(0), compare (J17l) . and to solve the equivalent equation 
of motion 



dv 



\^ 



'^ i=0 



(39) 



with i^(7o) = vq := cos{9q) and sgn(^(7o)j = <7i, := —ae- If Qvi^) has a zero of multiplicity two or more the 
solution can be solved by elementary functions. In general Qu{v) is a polynomial of order four with simple zeroes 
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(a) eQ : 



(b) eQ = 10 



FIG. 



Orbit configurations for the radial motion with a — 0.6, K — 1, V" — 0.4, and 1 < \/ -^ttvt < eQ. For a 



general description of colours and linestyles see the text. The spherical orbits on the blue (green) solid lines starting 

at a dot and approaching the light blue dashed (black dash dotted) lines are stable. In addition, the spherical orbits 

on the lower red solid line are stable but all others are unstable. Small plots on the top are enlarged details of the 

lower plot. 



only and can be solved with the following procedure that uses the Weierstrass elliptic function p (see [19|, |21|] . We 
transform the equation to the Weierstrass form by the substitution i/ = (— ^ — t^)^^ + ^6i where vq is an arbitrary 



zero of Q,^ and a, = — - — — 



(4-j)! djyC-i- 



ff (i^e)- This leads to 



4?^ - 992^ - ges 



(40) 



with 



982 



12 



2 1 
"2 - Tflias, 903 = 



216 



48 



010203 



16 



aoOg, 



(41) 



The initial conditions are ^(70) = Co '■= \{ ^°^^ + ^) and sgn ( ^(70) ] = (^i '■= ^^S^{o-z)^v The solution of eq. PO)) 
can now be expressed in terms of the Weierstrass elliptic p function. 



^(7) = p{l-le,in\ 982^993), 



(42) 
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K = 0.08 



K = 0.6 



K = 1.0 



FIG. 9: Orbit configurations for the radial motion with a = 0.6, P^ = 0.4, eQ = 0.3 and varying K. For a general 
description of colours and linestyles see the text. Note that at ^ = 0.0684 the plots change from case (ii) to (i). 




p2 =0.1 



p2 =0.195 



p2 = 0.2 



FIG. 10: Orbit configurations for the radial motion with a — 0.6, K = 0.8, eQ — 1.2 and varying V^. For a general 
description of colours and linestyles see the text. Note that at V^ ~ 0.1956 the plots change from case (in) to (iv). 



V^ = 0.5 



V^ = 0.627 



:0.63 
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FIG. 11: Orbit configurations for the radial motion with a — 0.6, K — 0.8, eQ — 0.9 and varying V^ . For a general 
description of colours and linestyles see the text. Note that at V^ ~ 0.6277 the plots change from case (i) to (ii). 



where 79411 is a constant such that p(7o — 70, in) = Co and sgn(p'(7o — 7e,in)) = fj- This is e.g. fulfilled by 704^ 



70 



-C/C 



?0 



dr 



V*^ 



-902^-903 



(with the principal branch of the square root). The solution for 9 is then given by 



0(7) — arccos 



03 



4p(7-70,m;5'02;.903)- -T 



ve 



(43) 



Let us now consider orbits which reach the poles. First, let 6'(7) be in the open interval (0, tt) for 7 e (71, 72) but on the 
endpoints 71^2 the orbit may reach the poles, 6'(7i_2) € {0, tt}. Then the solution on (71, 72), given by eq. (H51) . is also 
valid on the complete interval [71 , 72] because the right hand side of eq. (j331) is continuous on the whole closed interval 
with limits ^(71, 2) as 7 approaches 71,2- In general let 7^, i > 1, be the parameters with ^(7^) € {0,7r} and 7^ < 7^+1. 



Define 6',- 



l[7.- 



yA and solve the differential equation in each interval with the condition ^^(7^-1) = 0.i_i(7i. 



sgn ( ^(7i_i) j = — sgn ( —5-^(74-1) ) • The switch in sign of ^ in 7.; canonically identifies 9 on [0, vr]. 



B. r- motion 



The procedure to solve the equation of motion for the r, see P^ 

(^\ = n\f) - (ef2 + K)K{f) = R{f) 



(44) 



is analogous to that in the previous section. Again, the right hand side is a polynomial of fourth order. If it has a 
zero of multiplicity two or more the differential equation can be solved in terms of elementary functions. The general 
case can be solved with the substitution f = C3 (4^ ~ ^) + ^-R) where fn is a zero of R and Ci = (4^.a! %(i-i) {^R)- 
This leads to 



"(7) = 



C3 



4p(7-7r,in;5r2;5r3)- f 



rR 



(45) 



with gf2, 3f3 given as in (|4ip with a^ = q and gf^ — gg^. The parameter 7f,in only depends on the initial 
conditions, p(7o - 7r,in;ffr2;5r3) = l{nr^.+f) ^^'^ sgn(p'(7o - 7p,i„;5f2;5f3)) = -sgn(c3)CTp, e.g. 7p,in = 
7o + sgn(c3)CTf /_^° 



dT 



y/4T'-^-gf2'r-gf3 
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-motion 



The equation describing the (/)-niotion, see (|13p. can be rewritten as 



,(,).,„+ r^.,-riw,,, ,46, 

J^g A(r) J^^ sm-'6' 

where the right hand side is separated in a part which only depends on f and one that only depends on 0. We will 
now treat both integrals separately. 

The 8 -dependent integral 
Let us start with the integral 

Ueh)=r^dj. (47) 

J-yo Sill ^ 

If we insert the expression for 9 given by eq. (j43l) . which we write symbolically as 9 — 9(p{'j — "fg^in)), we get 



70 sin^ ^(p(7 - 7e,in)) ^o 



^^.(7)=/ ,'\^';.. ■'i. drf^ I i?0,(p(7-7.,in))rf7 (48) 



with R^{9(p{^ — 7e,in)) a rational function of p(7 ~ 79, in)- A partial fraction decomposition then yields 

l4>e (7) = r<^o + -. "'\ , + -, "'\ . rf7 (49) 

J70 p(7-7e,in) -Pei p(7 - 7ean) - Pe2 

with 

_ ePvQ-L l(eP±Z)a4 1 804 ± 03 - 1^903 /^^n 

The integral over each summand in eq. (|49l) can be expressed in terms of the Weierstrass Q = C('i5e2i3S3) ^-^i^ 
CT = ct(-; 5^2; flea) function (see [ll,[2l|) 

r / N ^ ^ , aei /, cr(7-«i) , ct(7 + 'Ui) 

^00 (7) =ae(7 - 70) + —rr^ [ In -7 r - In —, — - — r + 2(7 - 7o)C(wi) 

ae2 ^ g(7 ~ V2) , cr(7 + ^2) , „, .^, A ,^,^ 

'In— 7- In— ■ r +2(7-7o)C(w2) (51) 



p'{v2) V cr(7o - V2) cr(7o + 1^2) 

where Wi,f2 need to chosen such that p(ui + 79, in) = Pei, and p(w2 + 7e,in) = /302- 

T/ie f -dependent integral 
The procedure to solve the f-dependent integral 

r-an{f) 

is analogous to the previous section. We substitute f = f(p(7 — 7f,in)) of eq. (|45p and obtain a rational function as 
integrand, R^^ of p(7 — 7f,in), which we decompose into partial fractions 

I^. (7) = r R<pM^ - 7r-,in))d7 (52) 

70 P(7) - Pfl P(7) - Pr2 
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Here the constants are given by 



and ari 2 the coefRcients of the partial fraction decomposition. The solution has the the form of eq. ([5T|) with 
aei,2 = afi,2, fisi.2 = Pri,2 and wi,W2 such that p{vi + 7f,in) = /3ri and p(i;2 + 7r,in) = /3f2- Also the p, C, and cr 
functions refer to the parameters gr2i 9r3- 

D. t-motion 

The solution to the differential equation describing the t-motion, see (I14p . can be found in a similar way as for the 
(j) motion. An equivalent formulation of the differential equation is 

til) ^to+ r (""'+°y("") rf^ _ r -ari9)d^ , (55) 

where the right hand side is a again separated in an f and a 9 dependent part. We treat these integrals separately. 

The 9 -dependent integral 

We solve the integral 

Ite{l)-= f -ane)d^ 

with the same ansatz as in the foregoing section: First, we substitute 9 = 0(p(7— 7e,in)) of eq. p3|) and then decompose 
the integrand Rt^ which is a rational function of p(7 — 7e,in) in partial fractions, 

ri 

Iteil) = i?te(p(7 - 7e,i„))d7 

70 P(l) - Po (p(t) - Pef 

where 

ag — —aL — cl^Euq + a E + aePiyQ , (3$ = —03 

agi = --dai{2aEvQ - eP) , 002 = --rp:^'^^^^: ■ 

4 io 

The solution written in terms oi Q — C{'t 962t 983) and a — <y{-', ge2i 993) is given by 

cr(7 — v) a{'-f + v) 



(57) 



7 = «e 7 - 70 + ^ In -^ '-Z ~ In }'\ + 27 " 7o C ^^ 

p'('y) V cr(7o-w) cr(7o + w) 



a92 



C^92 ( 

727T - C(7 - «) + C(7o -v)- C(7 + w) + C(7 + w) - p(w)(7 - 7o) 



^ P'2(t;) 



2p' 
where v such that p(w + 7e,in) = P>9- 



Now the Integral 



p"{v) cr7-w cr(7 + w) , ^, \r( \\\ fr.Q\ 

77T(ln-7 T^l'^-? — — ^+2(7-7o)C(«)) ) (58) 



The r-dependent integral 



Jlo ^(^) 
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FIG. 12: Charged particle orbit in Kerr-Newman space-time. Here a B+ orbit in the northern hemisphere (A'" orbit) 

with a = 0.6, P = 0.47, Q = 0.16, e = 158.11, E = 1.98, L = -62.62, K = 33.33 is shown. Left: three-dimensional 

plot. The gray spheres correspond to the horizons. Right: projection on {x,z) plane. The dotted lines correspond to 

extremal 9 values. 



will be solved. The substitution r — r{p{j — 7f,in)) of eq. ([IS]) again leads to a rational function Rtf of p(7 — 7f,in), 
which becomes after a partial fraction decomposition 



itf 



R4'fipii-if:m))d-f 



70 

7 



af2 



p{l) - /3ri p(7) - /3f : 



-dj . 



(60) 



Here a^, /3ri 2 ^-re defined as in eq. (|54p . and a? 1.2 £^re the coefficients of the partial fraction decomposition. The 
solution has the the form of eq. (|5ip with agi2 — ari.2; f^ei 2 = /?ri 2 ^'ii'i ''^1)^2 such that p{vi + 7f,in) = /3fi and 
p{v2 + 7f,in) = I3f2- The p, C, and a functions here refer to 3^2 and 17^3. 



E. Examples 



The analytical solutions to the equations of motions are given by (H5)l , (P5|) , ([5T|l , and (155)1 , respectively with the 
appropriate constants ([54)1 . Here, we use these results to exemplify the orbit structure in Kerr-Newman space-time, 
see Figs, [m and [H 



F. Periastron shift and Lense-Thirring effect 



In General Relativity bound orbital motion can be much more complicated then the closed ellipses of Newtonian 
Gravity. However, in the weak field they are quite similar and the deviation can be characterized by a precession of 
the orbital ellipse, called the periastron shift, and the orbital plane, called the Lense-Thirring effect. They are caused 
by a mismatch between the periodicity of r{ip) {0{ip)) with 27r. These notion may be generalized to orbits in the 
strong field as demonstrated by Schmidt Q and Fujita and Hikida [lg| for Kerr space-time. An analogous treatment 
is also possible in Kerr-Newman space-times. 

For bound orbits, the radial and colatitudinal periods yl,. and Ag with respect to Mino time are defined by the 
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FIG. 13: Charged particle orbit in Kerr-Newman space-time. Here a F^ orbit in the northern hemisphere {N orbit) 
with a = 0.6, P = 0.47, Q = 0.16, e = 158.11, E = -11.83, L = -5.18, K = 33.33 is shown. Left: three-dimensional 
plot. The gray spheres correspond to the horizons. Middle: projection on (a;, z) plane. The dotted lines correspond 

to extremal values. Right: projection on (x, y) plane. Dotted lines correspond to the horizons. Note that 0(7) 
diverges at the horizons at some finite 70. In this plot we stopped at some 70 — A.y (A.^ small) and continued with 



70 ■ 



A^ on the other side of the horizon. 



smallest non-zero real value with r{\ + A^) = r(A) and ^(A + Ag) ~ 6{\) giving (see (|12l) and (TTTI) ) 

df . /■"""'- dO 



Ar = 2 



^p vW) 



Ab = 2 



VW)' 



(61) 



where rp is the periapsis and ra the apoapsis. To determine the periodicity with respect to if, we need to know the 
secular accumulation of (p with respect to the Mino time. This can be achieved by setting 



where ^{r,9) — $r(r) + ^e{0) is the right hand side of p^ . Here 



T 



mr,0))^:-- 



lim 



(A2-Ai)^cx) 2(A2 — Ai) J)^^ 



$(r, 6) dX 



(62) 



(63) 



is an infinite Mino time average and $osc(^) represent oscillatory deviations from this aver age. Using the averaged 
ip{X) = Ti^A, the periodicity of r{(f) {0{(p)) is then given by T^pAr {T^pAg). Accordingly, the periastron shift and the 
Lense-Thirring effect per revolution can be computed by 



Ap = T^ylr-27r, 



A 



LT 



T„yl« - 27r . 



(64) 



For neutral test particles and small a the periastron shift and the Lense-Thirring effect was calculated to first order 
in m. 



SUMMARY AND CONCLUSIONS 



In this paper we discussed the motion of charged particles in the gravitational field of Kerr-Newman space-times 
describing stationary rotating black holes with electric and magnetic charge. We demonstrated that it is sufficient to 
consider test-particles with electric charge only as an additional magnetic charge would only lead to reparametrization. 
After that we classified the orbits in radial and colatitudinal direction. For both a large variety of orbit configurations 
was identified as summarized in tables |T] and [Hi In particular, we also identified orbits crossing the horizons or r = 0. 
These configurations were then assigned to regions in the parameter space pictured in several figures. The boundaries 
of these regions are amongst others given by parameter combinations which represent orbits of constant r or 9, which 
were discussed in detail. For all orbit configurations analytical solutions to the equations of motion were presented in 
terms of elliptic functions dependent on the Mino time. 
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For the sake of completeness we considered here a black hole endowed with magnetic charge although such was not 
found until now. This has a big impact on the motion in direction in contrast to electric charge, which does not 
influence the colatitudinal motion at all. Not only deviates the motion from the symmetry to the equatorial plane for 
a nonvanishing magnetic charge but also additional types of orbits appear. For example, stable off-equatorial circular 
orbits outside the horizon do exist in this case, which are not possible else |l2l]. They are given by the intersection 
of a red solid line with _E^ < 1 as in Fig. [T]and a green solid line with £'^ < 1 as in Fig. [5l (E.g. a — 0.5, Q = 0.3, 
V = 0.6, eP « 2.1, K = 1, E Ki 0.9, and L sa 2.4 results in an orbit with constant r = 3 and cos9 w 0.8.) On the 
contrary, the magnetic charge does not influence the radial motion as is appears only in the combination P^ + Q^. 
Still, to our knowledge the discussion of the radial motion in this paper is the most complete so far for Kerr- Newman 
space-times. 

The analytical solution presented here are largely based on the 19th century mathematics of elliptic functions 
already used by Hagihara [16| to solve the geodesic equation in Schwarzschild space-time. However, a key ingredient 
here is the introduction of the Mino time which decouples the radial and colatitudinal equations of motion. We 
presented the results here in terms of Weierstrass elliptic functions which may be rewritten in terms of Jacobian 
elliptic functions. The advantage of our presentation is that one formula is valid for all orbit configurations. 
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